Viscosity and dilepton production of a chemically equilibrating quark-gluon plasma at 

finite baryon density 
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By considering the effect of shear viscosity we have investigated the evolution of a chemically 
equilibrating quark-gluon plasma at finite baryon density. Based on the evolution of the system 
we have performed a complete calculation for the dilepton production from the following processes: 
qq—*ll, qq^gll, Compton-like scattering (qg—*qll, qg—>qll), gluon fusion gg—*cc, annihilation qq— *cc 
as well as the multiple scattering of quarks. We have found that quark-antiquark annihilation, 
Compton-like scatterring, gluon fusion, and multiple scattering of quarks give important contribu- 
tions. Moreover, we have also found that the dilepton yield is an increasing function of the initial 
quark chemical potential, and the increase of the quark phase lifetime because of the viscosity also 
obviously raises the dilepton yield. 

PACS numbers: 12.38.Mh, 25.75.-q, 24.10.Nz 



I. INTRODUCTION 

The Relativistic Heavy-Ion Collider (RHIC) at the 
Brookhaven National Laboratory and the Large Hadron 
Collider (LHC) being built at CERN will provide the 
best opportunity to study the formation and evolution of 
quark-gluon plasma (QGP). Dileptons have large mean 
free path due to the small cross section for electromag- 
netic interaction in the plasma, therefore can provide an 
ideal probe for the detection and study of the plasma. 

Many authors [1-3], considering that the created QGP 
in collisions is a thermodynamic equilibrium system, have 
studied the dilepton production. Recently, the pho- 
ton and dilepton productions were studied based on the 
evolution model of chemically equilibrating QGP, estab- 
lished by Shuryak, Biro and co-workers [4-6], especially, 
the studies of these productions in the plasma at fi- 
nite baryon density were also performed [7,8]. However, 
most previous works were done by assuming the par- 
tonic plasma to be ideal, i.e., without any viscous effect. 
In principle, the viscous effects in fluid hydrodynamics 
should not be neglected in a realistic scenario since the 
dimension of the plasma is comparable to the mean free 
path of the partons. The viscous coefficient in the frame- 
work of hydrodynamics is composed of bulk and shear 
viscosity, while the bulk viscosity vanishes for a quark- 
gluon plasma [9]. In this work, we mainly discuss the 
effect of the shear viscosity, which have attracted many 
authors to investigate its influences on the formation and 
evolution of the QGP system. Authors of [9] have studied 
the viscous corrections to the hydrodynamic equations 
describing the evolution of the QGP at finite baryon den- 
sity, and investigated the effect of viscosity on chemical 
equilibration of the system. They have found that due 
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to the viscosity the lifetime of the plasma increases, the 
temperature evolution of the system becomes slow, and 
the chemical equilibration of the system becomes fast, 
therefore, the reaction rate will be heightened. How- 
ever, we should point out that in previous work many au- 
thors have regarded the viscous coefficients as adjustable 
parameters [9-12]. Indeed, they should be directly ob- 
tained from the thermodynamic quantities of the sys- 
tem. On the other hand, the viscous coefficients derived 
by Danielewicz and Gyulassy [13] based on QCD phe- 
nomenology for a baryon free-plasma, especially by Hou 
and Li [14] considering the Debye screening and damp- 
ing rate of gluons for a baryon-rich plasma using finite- 
temperature QCD, are so large that the temperature of 
the plasma would be abnormally heightened. 

In early calculations, one mainly considered the dilep- 
ton production from the process qq—>U. In recent years, 
possible sources of dileptons, such as qq—>ll annihila- 
tion, qg—>qll Compton-like scattering and qg—>qgll fu- 
sion, were investigated [15,16]. In addition, the contribu- 
tions of gluon fusion gg—>cc, quark-antiquark annihila- 
tion qq—>cc and multiple scattering of quarks to dileptons 
have also been studied [17]. 

In this work, starting from the shear viscous coeffi- 
cient given by relativistic kinetic theory for a massless 
QGP under relaxation time approximation, we first es- 
timate the mean free paths of partons in a chemically 
equilibrating QGP at finite baryon density, then com- 
bining with the parton energy densities, calculate the 
shear viscous coefficient of the QGP. Subsequently, based 
on our evolution model including the viscosity, we per- 
form a complete calculation for the dilepton production 
from processes: qq—>ll, qq^gll, Compton-like (qg^qll, 
qg—>qll), gluon fusion gg^cc, annihilation qq—>cc as well 
as multiple scattering of quarks to predict the contribu- 
tions of these reaction processes, and reveal the effect 
of the finite baryon density and viscous phenomena on 
dilepton production. 
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The rest of the paper is organized as follows: Sec. II 
describes the evolution of the dissipative QGP system. In 
Sec. Ill, we discuss the yields of dilcptons of the system. 
We give the results and discussions in section IV. Finally, 
in Sec. V, we give brief summary and conclusion. 



II. EVOLUTION OF THE DISSIPATIVE QGP 
SYSTEM 

In this work, we describe the distribution 
functions of partons with Jiittner distributions 
fq(q)= X q(q)/( eipT ^ q)/T + X q(q)) for q uarks (antiquarks) 
and f g —\g/(e p / T — X g ) for gluons, where fugacity \ t 
(<1) of the parton of type i is used to characterize the 
non-equilibrium of the system. Based on these distri- 
bution functions, we first derive the thermodynamic 
relations of the chemically equilibrating QGP system 
at finite baryon density. Expanding densities of quarks 
(antiquarks) 



n q(q) = 



9q(q) 

2ir 2 



9(9) 



/ 



p 2 dp 



(1) 



over quark chemical potential fj, q , we get the baryon den- 
sity of the system [18] 

n ^q = ^ T ^Ql\ - Qi x q) + ^ q T 2 (Q\K + Q\Xq) 
+T^ q {Q\X q -Q\X q ) + + (2) 

and the corresponding energy density 

e QG p = ^[T\Q\\ q + Q\X q ) + 3 N T 3 (Q 2 X q - Q 2 X q ) 



+3ii 2 q T 2 (Q\X q + Q\X q ) + Tfj? q {Q\X q - Q\X q ) 



where g q ^ and g g are degeneracy factors of quarks (anti- 
quarks) and gluons respectively. Since the convergence of 
the following integral factors appearing in the expansion 
above 



rr, 



Z n dZ 



Z n dZ 



Z n dZ 



(e z + x q y 
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is very rapid, it is easy to calculate these integral numer- 
ically [18]. 

We consider the reactions leading to chemical equilib- 
rium: gg^ggg and gg^qq. Assuming that elastic par- 
ton scatterings are sufficiently rapid to maintain local 
thermal equilibrium, the evolutions of gluon and quark 
densities can be given by the master equations, respec- 
tively. We first extend the master equations to include 
the viscosity as done in [9]. Similarly, the evolution of 



baryon density can be described by a corrected conser- 
vation equation of baryon number including a viscous 
term. In addition, due to viscosity, a viscous term would 
be contained in the conservation equation of the energy- 
momentum, too. Combining the master equations to- 
gether with the equation of baryon number conservation 
and equation of energy-momentum conservation includ- 
ing viscous corrections, for longitudinal scaling expansion 
of the system, one can get a set of coupled relaxation 
equations (CRE) describing evolutions of the tempera- 
ture T, quark chemical potential fi q , and fugacities X q 
for quarks and X g for gluons on the basis of the thermo- 
dynamic relations of the chemically equilibrating QGP 
system at finite baryon density [9,18] 
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where n q i q \ is the value of n„c«i at A„r^i = 1, n„ 
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l q(q) 
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n q /(g q /2Tr 2 ), n° = n g /(g g /2n 2 ), £(3)=1.20206, and r? the 
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shear viscous coefficient. The gluon and quark produc- 
tion rates R3/T and R2/T are respectively given by[6,18- 
20] 



3ai A g L 9g 2 T 4 /2 



R 2 /T = 



24tt G 2 



iV^ln^), 



(10) 



7T J 



(11) 



where M|, is the Debye screening mass, g 2 — Aira s 



and I(\ g , \ q ,T, /ig) is the function of X g , X q , T, fi q , as 



used in [5-6]. We here take the quark flavor A/=2.5 [18- 
20]. Solving the set of evolution equations (5) — (8) under 
given initial values obtained from Hijing model, we can 
get the evolutions of temperature T, quark chemical po- 
tential ii q and fugacities \ for quarks and X g for gluons. 

To discuss the effects of shear viscous coefficient, we 
have quoted two different expressions of it: 771 and 772 are 
taken from [9,13], respectively. 



m = m 



e QGP 



T 

112 = ~ 



a v jn g + n q gU q + n g 



(12) 



(13) 



where 770 is treated as a constant [9] , and a v the transport 
cross section [13]. 

Now, we discuss the calculation of the viscous coeffi- 
cient rj in our work. According to [13,14,21], the shear 
viscous coefficient using the relativistic kinetic theory for 
a massless QGP in the relaxation time approximation is 
written as: 



(14) 



where Aj is the mean free path of particle of type i in 
QGP, which in a chemically equilibrating QGP for gluon 
is given by 

4 1 Ml{Ml+$T 2 /2) 



^ 9 9n n 27ra s 2 



9T 2 /2 



and for quark by 



A„ 



1 M 2 D {M 2 D +9T 2 j2) 



An q 2-KOt s 2 



9T 2 /2 



(16) 



where M|> is the Debye screening mass and given by (11). 
Thus, we can directly calculate the viscous coefficients r/ g , 
i] g and their total 77 from the thermodynamic quantities 
of the QGP system. 



III. 



DILEPTON PRODUCTION 



Based on the evolution of the QGP system, we first 
consider dilcpton production from quark annihilation 
qq—>gll and Compton scatterings qg^qll and qg^qll. 
Their production rates can be calculated by [16] 



E 



dR 

d 3 p 



1 



/ 



d 3 pi d 3 p 2 d 3 p 3 
2(2tt) 8 J 2E X 2E 2 2E 3 

h(E 3 ))-5 4 (P 1 +P 2 -P 



f 1 (E 1 )f2(E 2 )(l± 



(17) 



where f(E) is the Jiittner distribution function of par- 
tons, ^2 I M 1 2 the square of the matrix element for re- 
action processes summed over spins, colors and flavors. 
The plus sign is for the annihilation process and the mi- 
nus for the two Compton processes. According to [22], 
above equation can be rewriten as 



dR 

d 2 M 



100 a 2 a s 
~2H tt 5 M 



dsdt 



dE x 



u 2 +t 2 + 2sM 2 
ut 
dE 2 



X, 



e (£i-^)/T + Xq e (Ei+n q )/T + Xg 



A„ 



1 6(P(E 1 ,E 2 )) 

for the annihilation process, and 



(18) 
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= 20 
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"AgAg 



e (E lT „ q )/T + Xq e ( El )/T _ X( 

\ 0(P(E 1 ,E 2 )) 

ei E 1+ E 2 -E T „ q )/T + \ q \p(E ll E2)) 1 / 2 



(19) 



for Compton scatterings, where P{E\,E2) = —{tE\ + 
(s + t)E 2 ) 2 + 2Es((s + t)E 2 -tEk)- s 2 E 2 + sH + st 2 ,9 is 
the step function, a the fine-structure constant, and a s 
the running coupling constant. The minus sign is for the 
Compton process qg—>qll and the plus for qg—>qll. The 
letters s, t and u are the Mandelstam variables. And the 
integration are performed over — (s — M 2 ) + k 2 < t < —k 2 
and M 2 + 2k 2 < s < 00 [16]. The cutoff k 2 is replaced 
by the thermal quark mass 2m 2 , [22]. For a chemically 
(15) equilibrating QGP system at finite baryon density m 2 is 



given by 



4a., 



m 2 = ^T 2 [2(G}A 3 



>i\) + (^) 



A + 1 J 



(20) 



where integral factors G\ and Q\ have been given above. 

Obviously, above calculations give up the infrared con- 
tribution because of introducing the infrared cutoff k 2 . 
Authors of [22] have discussed the infrared contribution 
to photon production. Following their calculation, in this 
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work, we have given an assessment of the contribution 
from the infrared part to the dilcpton produnction. The 
dilepton production rate with total energy E and total 
momentum p can be calculated by [22,23] 



dR = a 1 1 i m n " 

dEd 3 p 12tt 4 P 2 e E / T - 1 R ^ 



(21) 



where Yl^ is the retarded photon self-energy. The in- 
frared divergence mentioned above is caused by propaga- 
tion of soft momenta. To cure the problem it is necessary 
to dress one of the quark propagators (as done in Fig. 3 
of [22]). Then the retarded photon self-energy can be 
written as [22] 



3 rj 



d 3 k 
(2nf 



Tr[S*(k)^S(p-k) 7 »], (22) 



where S*(k) is the dressed propagator for a quark with 
four-momentum k, and S(q) is the bara propagator for 
a quark with four-momentum q = p — k. Following [22], 
after performing some calculations and applying the el- 
egant method developed by Braaten, Pisarski and Yuan 
[23] , one can finally obtain a simplified expression of the 
imaginary part of the retarded photon self-energy 

Imn ^ = ^(^ /T -1) 

rk c pk 

/ dk duf(w)f(E -w){k-oj)(3+(w,k),(23) 

JO J-k 

where / is the Jiittner distribution function of partons, 
and 



P+(uj,k) = 



\m\(k 



{k(co - k) - m 2 q [Q Q {z) - Qi(z)}} 2 + [^2(1 - z)] 2 ' 

(24) 

where z = u/k [24], Q and Qi are the Legendre func- 
tions of the second kind. After performing integral over 
the dilepton energy (E > M) one can obtain the dilepton 
production rate over square invariant mass of dilcptons 
dR/dM 2 . 

Similar to the preceding treatment, the production rate 
of quark-antiquark annihilation qq—>ll can be given by 
[25] 



dR 

dM 2 



24ir 4 



M 2 a lT (M 2 ) 



OO 



oc 



< / d Pl f q { Pl ) / dp 2 f q (p2), (25) 

10 JM 2 /4 Pl 

where a^M 2 ) = ^47ra 2 /3M 2 is the quark annihilation 
cross section. 

Aurenche, Gelis, and their co-workers have studied the 
dilepton production from bremsstrahlung and off-shell 



annihilation where the quark undergoes multiple scatter- 
ing in the medium as shown in Fig.l of [26]. Their cal- 
culation includes the Landau-Pomeranchuk-Migdal effect 
and is concluded that this contribution is important. Ac- 
cording to these authors' approach, the imaginary part 
of the retarded current-current correlator H R ^ (P) in 
Eq.(21) can be computed by 



5 

6ir 



+oc C d 2 a 

dqo[f(k )-f(q )}xRe 1 11 

■oo 



la + „ f („ \ 



1 



2(qoh) 2 



\/\qoko\ V 



(2^) 2 
(26) 



with k = qo + E, f is the Jiittner distribution function 
of partons again, and the dimcnsionless functions f(q±) 
and g(q±) respectively obey the integral equations [26] 



2g ± = iSEf(q ± ) + -g 2 T 



I 



rC(i±) 



(27T) 2 

x[f(q±)-f(q±+l±)\ 



(27) 



and 



= iSEg( qi _) + g 2 C F T J |^jC(*±) 



(2n) 2 

x[9(q±)-g(q± + i±)}- (28) 

Using the method described in [26], we recast (27) and 
(28) as differential equations and solve them using a sim- 
ple algorithm, finally, can get Re J d 2 q±/ (2Tr) 2 q± ■ f(q±) 
and Refd 2 q±/(2ir) 2 g(q±) and the corresponding dilep- 
ton production rate of multiple scattering process. 

For the QGP system, produced in collisions at RHIC 
energies, with very high initial temperature («0.57 GeV) 
[6,18], thermal charmed quark production and its contri- 
bution to lepton pairs should be contained, especially, 
those from the gluon fusion gg—>cc and quark-antiquark 
annihilation qq^cc. Similar to the calculation for qq—>ll, 
replacing the cross section tT w -(M 2 ) appears in the expres- 
sion (25) with those of the reactions qq—>cc and gg-^cc 
in leading order QCD, we can compute the yields of 
charm pairs in the QGP. Almost all of the produced 
thermal charmed quarks would eventually hadronize to 
D-mesons [17]. Considering that the D-meson decays 
to leptons with a 17% branching ratio for charged D- 
mesons [17,27,28], finally one can obtain the contribution 
of charmed quarks from reactions gg^cc and qq-^cc to 
lepton pairs. 

We integrate these production rates over the space- 
time volume of the reaction. According to Bjorkcn's 
model, the volume element is d A x = d 2 XTdyrdT, where 
t is the evolution time of the system and y the rapidity 
of the fluid element. We consider Au 197 + Au 197 central 
collisions, so the integration over transverse coordinates 
just yields a factor of d 2 xr = kR 2 a , where Ra is the nu- 
clear radius. Finally we obtain the dilepton spectra of 
the system 



dN 
dydM 2 



ttR 2 , 



rdr 



dR 
dM 2 ' 



(29) 
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FIG. 1: (Color online) The evolutions of the temperature T with 
proper time r for viscous coefficients rji , r/2 , and r) at initial quark 
chemical potentials ^,,0=0.000, 0.284, 0.568, and 0.852 GeV. The 
rji is the viscous coefficients for parameters 770=0.0, 0.4, and 0.8 as 
given in [9], and the corresponding values are denoted by the solid, 
dash and dot curves, respectively. 



IV. CALCULATED RESULTS AND 
DISCUSSION 

In this work, we focus on discussing Au 197 + Au 197 
central collisions at the RHIC energies. In order to com- 
pare with [9,29,30], we take initial values of the system: 
t =0.70 fm, T o =0.570 GeV, A g0 =0.08 and A g0 =0.02 from 
the Hijing model calculation. We have investigated the 
effect of viscosity on the evolution of the system us- 
ing those expressions of the shear viscosity as shown in 
(12) — (16). To understand the effect of the baryon den- 
sity on the dilepton production, we have solved the CRE 
for initial quark chemical potentials /u g o=0.000, 0.284, 
0.568, and 0.852 GeV for given viscous coefficients, and 
obtained the evolutions of the temperature T, quark 
chemical potential fi q and fugacities A g and \ q of the 
system. In Fig. I, we have shown the evolutions of the 
temperature T with proper time r for the various viscous 
cofficients from the expressions in [9,13]. Following [9], 
we take the parameters 770=0. 0, 0.4, and 0.8. The corre- 
sponding curves are indicated by the solid, dash and dot 
lines, respectively. We can see that the viscosity leads to 
the increase of the temperature of the QGP system. We 
also note that the temperature is unreasonably height- 
ened for viscous coefficient 772. For 771, the evolution of 
temperature seems to be reasonable, however, the viscos- 
ity rji is only obtained through adjusting the parameter 
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FIG. 2: (Color online) The viscous coefficient r\ as a function of 
the temperature. Black, red, green and black denote, in turn, the 
calculated viscous coefficients r\ for initial quark chemical poten- 
tials ^ g o=0.000, 0.284, 0.568, and 0.852 GeV at the same initial 
conditions as given in Fig.l. 



t]q. In this work, we have calculated the viscous coef- 
ficients 77 by the thermodynamic quantities of the QGP 
system using (14) — (16). From Fig.l, we note that the 
calculated temperature distribution is reasonable. In ad- 
dition, the temperature is also a increasing function of the 
initial quark chemical potential. Fig. 2 shows the value 
of 77 as a function of the initial temperature, where the 
black, red, green and blue curves denote, in turn, the cal- 
culated 77 for initial quark chemical potentials /^ g o=0-000, 
0.284, 0.568, and 0.852 GeV at initial values mentioned 
above. From Fig. 2, one can see that 77 increases with 
increasing the temperature T. 

The estimated evolution paths of the system in the 
phase diagram have been shown in Fig. 3, where black, 
red, green and blue curves are, in turn, the calcu- 
lated evolution paths for initial quark chemical poten- 
tials /j 90 =0.000, 0.284, 0.568, and 0.852 GeV. The solid 
line is the phase boundary between the quark phase and 
hadronic phase. The lines with open circles denote the 
evolution of system without viscous effect, while the lines 
with solid circles are from the system with viscosity. The 
time interval between the two circles is 0.3 fm. The cor- 
responding equilibration rates of gluons and quarks, A g 
and X q , are shown in Fig. 4. The solid lines are for the 
cases of viscosity, and the short dash lines denote ideal 
cases, where the black, red, green and blue curves are, 
in turn, the calculated values for initial quark chemical 
potentials /j g0 =0.000, 0.284, 0.568, and 0.852 GeV. From 
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FIG. 3: (Color online) The calculated evolution paths of the sys- 
tem in the phase diagram for the initial values as given in Fig.l, 
where black, red, green and blue curves are, in turn, the calculated 
evolution paths for initial quark chemical potentials /x 9 o=0.000, 
0.284, 0.568, and 0.852 GeV. The lines with open circles denote 
the evolution of system without viscous effect while the lines with 
solid circles are the evolution of the system with viscosity. The 
time interval between the two circles is 0.3 fm. 
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FIG. 4: (Color online) The calculated equilibration rates at the 
same initial conditions as given in Fig.l. The solid lines are for 
the cases of viscosity while short dash lines for the ideal cases, 
where the black, red, green and blue lines are the calculated values 
for initial quark chemical potentials /^ q o=0.000, 0.284, 0.568, and 
0.852 GeV. 
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FIG. 5: (Color online) The calculated dilepton spectra from quark 
annihilation processes qq-^ll and qq— >gZZ. The signs of the curves 
are the same as those in Fig. 4. 

Figs. 3 and 4 we see that the evolution of the system be- 
comes slower due to viscosity, whereas the equilibration 
rate of the plasma becomes faster compared to the one 
in the ideal case. And also the effect of the initial quark 
chemical potential on the evolution is in accordance with 
previous conclusion [18]. 

Based on the evolution of the system described in 
Figs. 3 and 4, we have calculated dilepton spectra for pro- 
cesses qq^ll, qq—*gll, qg^-qll, qg—>qll, qq^cc, gg^cc, 
and multiple scattering at the four initial quark chemical 
potentials as given above. The calculated dilepton spec- 
tra from quark annihilation processes qq—>ll and qq^gll 
are, in turn, shown in panels (a) and (b) of Fig. 5. The 
solid lines are the spectra of the system with the effect of 
viscosity while dash lines indicate those from ideal sys- 
tem, where the black, red, green and blue lines denote 
the calculated spectra for initial quark chemical poten- 
tials ^ 9 o=0.000, 0.284, 0.568, and 0.852 GeV. One can 
see that the calculated spectrum goes up with increas- 
ing the initial quark chemical potential. The law is valid 
for processes qg^qll and multiple scattering as shown 
in panels (a) and (c) in Fig. 6. We know that the quark 
density goes up with the increase of quark chemical po- 
tential while the density of antiquark goes down, which 
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FIG. 6: (Color online) The calculated dilepton spectra from pro- 
cesses qg^qll, qg^qll and multiple scattering. The signs of the 
curves are the same as those in Figs. 4 and 5. 



will lead to a suppression of the production from process 
qg^qll, as shown in panel (b) of Fig. 6. And also we can 
see from Figs. 5 and 6 that the dilepton production is a 
increasing function of viscosity, which would be mainly 
attribute to the increase of the quark phase life-time due 
to the viscosity. 

The productions of soft dileptons, which are connected 
with the infrared contribution, have been computed fol- 
lowing the method represented in [23] for initial quark 
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FIG. 7: (Color online) The calculated spectra of all processes for 
the initial quark chemical potential /i 9 o=0.284 GeV. Curves 1—9 
are, in turn, the calculated spectra for qq— *cc, gg^cc, qg^qll, soft 
dileptons, multiple scattering, qg^qll, qq^gll, qq^ll, and their 
total. 
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FIG. 8: (Color online) The calculated total dilepton spectra of the 
system. The signs of the curves are the same as those in above 
figures. 
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chemical potentials m 9 o=0.000, 0.284, 0.568, and 0.852 
GeV. In Fig. 7 we have shown the results from all pro- 
cesses and their total for the initial quark chemical po- 
tential /i g o=0.284 GeV. Curves 1—9 represent, in turn, 
the calculated spectra for qq^cc, gg^cc, qg—>qll, soft 
dilcptons, multiple scattering, qg^qll, qq^gll, qq-^>ll, 
and their total. From Fig. 7, one can see that the spec- 
tra from the quark-antiquark annihilations qq—>ll and 
qq—>gll dominate. The infrared contribution is as im- 
portant as that of reaction qq^gll and even higher than 
the later one in the range of small invariant mass. The 
contributions from Compton-likc scattering qg—^qll, mul- 
tiple scattering, and annihilation gg—>cc can not also be 
neglected. 

We have also given the total yields of all processes 
of the system for initial conditions mentioned above, as 
shown in Fig. 8. The black, red, green and blue curves 
represent the total yields for ^0=0-000, 0.284, 0.568, and 
0.852 GeV, respectively. To understand the effect of vis- 
cosity on the dilcptons production, we also give the yields 
for ideal QGP system, which are denoted by dash lines. 
It shows clearly that the dilcpton yield of the system 
goes up with increasing initial quark chemical potential. 
However, previous authors have found that dilcptons pro- 
duced in a thermodynamic equilibrium QGP system are 
suppressed with increasing initial quark chemical poten- 
tial [1] . In this work, since that both the quark chemical 
potential and the temperature of the system are functions 
of time, compared with the baryon-free QGP it necessar- 
ily takes a long time for value {n q ,T) of the system to 
reach a certain point of the phase boundary to make the 
phase transition. Furthermore, in the calculation we have 
found that with increasing the initial quark chemical po- 
tential the production rate of gluons goes up, and thus 
their equilibration rate goes down, leading to the little 
energy consumption of the system, i.e., slow cooling of 
the system. Since gluons are much more than quarks 
in the system, with increasing the initial quark chemical 
potential the cooling of the system further slows down. 
These cause the quark phase life-time to further increase, 
as seen in Fig. 3. These effects will heighten the dilepton 
yield and compensate the dilepton suppression, leading 
the spectrum of the system to be an increasing function of 
the initial quark chemical potential. On the other hand, 
as seen in Fig. 3, due to the viscosity the evolution of the 
system becomes even slower, so that the dilepton yield 
will be heightened, as seen in Fig. 8. From Fig. 8 one can 
note that the dilcpton yields arc remarkably heightened 
due to the effect of the viscosity of the QGP system. 



V. SUMMARY AND CONCLUSION 

In this work, taking into account reactions gg^ggg 
and gg^qq leading to the chemical equilibrium of the 
QGP system, and conservations of energy-momentum 
and baryon number, as well as viscosity of the QGP sys- 
tem, we have derived a set of coupled CRE of the chem- 
ically equilibrating QGP system with viscosity at finite 
baryon density, produced from Au 197 + Au 197 central col- 
lisions at RfflC energies, which describes the space-time 
evolution of the system. Then, we have solved the CRE, 
and directly obtained the viscous coefficients from the 
thermodynamic quantities of the QGP system. We note 
that the calculated results of the viscous coefficients are 
reasonable. Subsequently, based on the evolution of the 
QGP system we have computed the dilepton spectra of 
the QGP system, we have found that the spectra is dom- 
inated by the quark-antiquark annihilation qq—>U and 
qq^gll, following by the multiple scattering of quarks, 
compton-like scattering qg^qll and annihilation gg^cc. 
While we have also calculated the infrared contribution 
and found it very important at the range of small in- 
variant mass of dileptons. Furthermore, we note that 
the increase of the dilepton yield with increasing the ini- 
tial quark chemical potential can compensate the dilep- 
ton suppression, thus eventually leading to the dilepton 
spectrum to be an increasing function of the initial quark 
chemical potential. Especially, we have found that the 
dilcpton yield of the system is obviously enhanced due to 
the viscous effect because this effect makes the evolution 
of the system slow down and thus the life-time of the 
QGP system increases. 
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